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1 Introduction 



The introduction of supersymmetry into classical mechanics has been motivated in part by the 
study of supersymmetric quantum mechanics, e.g. in the influential paper by Witten not 
the least with the aim to find a classical "background" for these models. 

That this can be achieved has been long known, namely that classical supersymmetric mechanics 
can be consistently constructed by replacing all dynamical quantities of the theory in question 
by even and odd elements of a Grassmann algebra B according to their respective bosonic or 
fermionic nature. This Grassmann-valued mechanics has been studied first by Berezin and 
Marinov M and Casalbuoni [[|. A more recent account devoted to one- dimensional systems 
was given by Junker and Matthiesen in M. For a general review of Grassmannian geometry 
the reader may turn to the book by de Witt ||. 

An interesting but somewhat different approach has been taken in a series of papers by Gozzi 
et. al. ||. Using path integral methods in classical mechanics he and coauthors construct 
an N = 2 supersymmetric theory for any given Hamiltonian system, where fermionic fields 
correspond to sections of various cotangent and tangent bundles of the configuration manifold. 
It should be noted, however, that the resulting model differs from the one analyzed in this 
paper, so that solutions or symmetries cannot be directly compared. 

While there is thus a considerable amount of literature on supersymmetric mechanics and 
interest in the topic has grown recently, there is a surprisingly low number of specific models 
that have been investigated and solved. Furthermore, attention is usually restricted to the 
case of one spatial dimension. An interesting exception on two-dimensional Liouville systems 
is given in @. 

We try in this paper to address both of these shortcomings, first by establishing some general 
properties of multidimensional supersymmetric mechanics, especially for centrally symmetric 
potentials and then by analyzing and solving a particular classical model, the Coulomb problem. 
The l/|a;|-potential is one of the oldest and most prominent potentials that has been studied 
in physics. This is not only since it properly describes both the gravitational and the electrical 
interaction of massive or charged point particles, respectively but also because it has some very 
attractive features. One example is that by a theorem of Bertrand || the attractive version of 
the potential is the only centrally symmetric potential next to the harmonic oscillator such that 
every admissable trajectory is closed. Furthermore one can derive an explicit time-dependent 
solution for every given set of initial data. This is a very useful property because it implies, 
if perhaps unexpectedly at first sight, that the supersymmetric version of the theory can be 
explicitly solved, too. 

This paper is based on the experience gained from an earlier investigation into one-dimensional 
supersymmetric problems. Both in |J and [|nj, where a different concept of complex conjuga- 
tion for Grassmann variables was used, one of the ideas of tackling the supersymmetric problem 
was the assumption of a finitely generated Grassmann algebra B. This led to a layer-by-layer 
structure of the theory which could be used to solve one-dimensional models by subsequently 
solving higher layers based on the solutions of all lower ones. We shall take a similar approach 
in this paper. 

The starting point for our investigation will be the statement of the Lagrangian in section 
0. It is obtained ultimately from the familiar N = 2 supersymmetric (1 + l)-dimensional field 
theory with Yukawa interaction via dimensional reduction which leads to the description of one- 
dimensional supersymmetric point particle mechanics analyzed in || . The resulting Lagrangian 
can then be generalized to more than one dimension. Notice that supersymmetry requires the 
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bosonic potential to be non-negative. This means that our attention has to be restricted to the 
case of repulsive interaction. 

We go on to derive the complete set of 3n equations of motion and show that as in the one- 
dimensional model there are two independent supersymmetries as well as a purely fermionic 
invariance transformation which can be seen as an internal "rotation" of the fermionic variables. 
We proceed by giving one explicit solution to the fermionic equations of motion utilizing the 
supercharges - except in the one-dimensional case this can not be the complete solution, though. 
We then restrict our attention to a special class of models, namely those with a centrally 
symmetric scalar potential. After restating the equations of motion we derive in section [I] the 
explicit form of the supersymmetric angular momentum and show that it is a supersymmetric 
invariant. 

The symmetries of our model can be understood best from its Poisson algebra which we therefore 
describe completely in section f|. 

In section |5| we write down the complete set of equations governing the supersymmetric Coulomb 
problem. Even in three dimensions these are still nine coupled first and second order differential 
equations and it seems by no means obvious that there will be an explicit analytical solution 
for every possible set of initial data. 

We use section |6| to show that such an explicit solution exists. This is achieved by utilizing a 
method originally employed in |L(| and then applied further in |J. The key idea is to decompose 
all dynamical quantities (and therefore every equation) into terms containing the same number 
of Grassmann generators. The model then naturally obtains a layer structure which can be 
analyzed from the bottom layer - the classical problem - upwards. In particular, solutions 
to higher layers may be found by using all those already derived for the lower layers. We 
will carry out this programme here up to second order in the Grassmann algebra - since our 
results indicate a close resemblance of the multi-dimensional to the one-dimensional case we 
can suspect that higher order solutions exhibit many of the characteristics described already 
for the one- dimensional case. We will give a short summary of these in our conclusions. 
As mentioned we begin our analysis in section |6.1| with the bottom layer of the system which 
is identical to the classical, i.e. non-supersymmetrized problem. Although the solutions to 
the problem are well-known we explain briefly how explicit time-dependent expressions can be 
found - we do this because the textbook treatments of the problem usually only derive the 
possible orbits but do not specify the actual functions of time that describe it. 
The fermionic equations of motion comprise the next layer and we study their solutions in 
detail in subsection |6.2j . As in the one- dimensional case we find that fermionic motion is not 
only compact but restricted to a spherical space. Furthermore, while there is a plane of motion 
for the bosonic variables (of lowest order) the same is not true for the fermionic variables or, 
in other words, fermionic motion is non-trivial in all three spatial directions. 
Finally, we solve the top layer of our system, consisting of the bosonic quantities of second order 
in the generators. Though the particular solutions found in the one-dimensional case cannot 
be applied directly to the Coulomb problem we can carry over the most important aspect to 
the three-dimensional case, namely that solutions to the homogeneous equations of motion of 
the top layer system are variations of the bottom layer system with the free parameters of that 
motion. In particular, we show that every top layer solution to the homogeneous equations of 
the Coulomb problem can be explained as a linear combination of the variations with initial 
time, energy, eccentricity, orientation of angular momentum and orientation of the hyperbola 
in the plane of motion. We conclude our analysis by a full statement of the explicit solutions 
of the inhomgeneous equations including the boson-fermion-interaction term. 



2 



As is well-known, the classical Coulomb problem exhibits a hidden 0(3, l)-symmetry as a 
consequence of which there is an extra vector- valued conserved quantity, the Runge-Lenz- vector. 
It comes as an unexpected surprise that this symmetry is broken in the supersymmetric version 
of the problem. We demonstrate this, i.e. that there can be no supersymmetric version of the 
Runge-Lenz vector, in the appendix. 



2 Multidimensional supersymmetric mechanics 

We assume that supersymmetric mechanical models can be described by the following La- 
grangian: 

c = \ & '* ~ \ U ^' U ^ + \^+'^+ + + #+VC/(x)V-, (1) 

which is a generalization of the familiar supersymmetric Lagrangian in one-dimensional me- 
chanics ultimately derived from (1+ 1)- dimensional field theory. The dynamical quantities 
x, tp+ and i/>_ are n-component vectors, e.g. x = (x±, . . . ,x n ), x is Grassmann-even and 
and are Grassmann-odd. 

In addition, U is taken to be a vector- valued n-component Grassmann-even function of x, 
so that VU is a n x n-tensor with elements diUj(xk)- We will usually take this tensor to 
be symmetric which will always be the case if U is derived from a scalar potential term as 
U(x) = VW(x). 

xy denotes the standard Euclidean inner product, extended if necessary in the obvious way: 
xMy, M being an n x n-tensor, stands e.g. for the component expression x l M i jy^ . 

The equations of motion can now be read off from the Lagrangian: 

x = -VU(x)U(x)+iV(iJ> + VU(x)ip-) (2) 
^+ = -VU(x)ij>- (3) 
^_ = VU(x)i/> + . (4) 

Every equation is vector-valued and has n components. 

Invariance with respect to time-translation means that there is a conserved Hamiltonian which 
can be calculated to be 

H = -x-x + hj(x)-U{x) - iil) + VU{x)ij)-. (5) 



In addition there are two independent supersymmetry transformations, namely 

5x = ieif>+, (5-0-1- = —ex, <5i/>_ = —eU(x) 1 
5x = ieij)-, 5tp + = eU(x), 5tp- = —ex, 



(6) 



where e is an infinitesimal scalar Grassmann-odd parameter. The corresponding conserved 
supercharges are given by the expressions 

Q = x-il>+ + U(x)-i>_ (7) 

Q = x-i/f--U(x)-ij> + . (8) 
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Finally, an internal transformation in the fermionic part of the Grassmann algebra, given by 



5ij) + = Trf-, Sip- = -rr^+, 

where rj denotes an infinitesimal scalar Grassmann-even parameter, also leaves the Lagrangian 
invariant. The resulting charge is simply 

R = iiJ> + .il>_. (9) 

One of the results from ||, that can be carried over from the one-dimensional case almost with- 
out alteration, is that solutions to the fermionic equations can be found using the supercharges. 
In one dimension fl7|) and (§) can be formally inverted to yield: 

V>+ = ^(Qx-QU) (10) 

</>- = ^(QU + Qx), (11) 

where E is the constant Grassmann-valued energy. This is still a solution in the n-dimensional 
case, though special care has to be taken if the real part Eq is zero. 

Unlike the situation in one dimension, however, (110) and ( |TT| ) give only part of the fermionic 
solution. In the mult i- dimensional case there is still a 2(n — l)-dimensional solution space that 
needs to be determined. 



3 Centrally symmetric potentials and angular momen- 
tum 

So far these formulae have been straightforward generalizations of formulae from the one- 
dimensional case. However, in more than one dimension there can be additional conserved 
quantities if there are extra symmetries in the system. We will now assume that the scalar 
bosonic potential 

V = l -U{x)-U{x) 

is central, i.e. a function of \x\ = (x-x)^ only. This can be easily achieved by demanding that 
the scalar potential W(x) is a function of \x\, for then U(x) will be of the form 

U(x) = W\\x\)x, 

with x denoting the unit vector £c/|cc|. 

The equations of motion can now be written as follows: 



x = -W'W"x + i[ Li- (V>+(V>--£) + (</>+•£)</'-) (12) 
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where the argument of W is \x\. Note that the purely bosonic part of the problem is described 
just by the first term on the right hand side of (|I!?D. 

For convenience we will now restrict ourselves to the three-dimensional case relevant for the 
Coulomb problem, although our discussion should generalize easily to more than three dimen- 
sions. 

The particular form of the potential function W means that there is an extra symmetry in the 
system, given by rotating both the bosonic and the fermionic variables by the same amount 
around the same axis a. Infinitesimally, we can write 

5 a x = e[a,x}; S a ij} + = e[a, V>+]; 5 a rj>- = e[a, t/>_], 

where [■, ■] denotes the three-dimensional vector product, i.e. [x, y] 1 = e^uXjUk- 
This leads to the conserved angular momentum: 

L = [x, x] - ~ V+] - \ [V>-, 1>-] • (15) 

Note that the second and third term of L are non-zero since if>+ and are Grassmann-odd. 
It is an interesting observation that the strict condition for the existence of angular momentum 
is that W is central. By choosing U = f(\x\) Ox; O G 0(3) we end up with a rotationally 
invariant bosonic potential V but angular momentum is not conserved. 

Finally we point out that the following two expressions are also time-independent: 

which can be said to denote the product of the oriented "areas" and "volumes" of and 
respectively. They do not, however, constitute independently conserved quantities as the first 
one is proportional to R 2 and the second one to R 3 . 

Apart from invariance under time translation, another important aspect is that of invariance 
under supersymmetry. It can be shown directly that angular momentum and energy are super- 
symmetric invariants whereas e.g. the extra charge R is not: 

5H = 0, 5L = 0, 5Q = -2eH, 6Q = 0, 5R = -xeQ 

5H = 0, 5L = 0, 5Q = 0, 5Q = -2eH, 5R = ieQ. ^ ' 

Here 5 and S are the two supersymmetry transformations introduced in equation ([]) and e is a 
Grassmann-odd infinitesimal parameter. 



4 Poisson- Algebra 

For a better understanding of the symmetries of our mechanical system we will now derive the 
Poisson-algebra of the problem. Therefore, we must write down the definition of the Poisson 
brackets for Grassmann-valued quantities. To do this we have to introduce the canonical 
momenta p and 7r by 

dL dL 3L 

dx' + dip+ dtp- 
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For the specific choice of our model one finds immediately that 



p = x, n + = --if> + , 7r_ = --V>_. (17) 

The fact that the canonical momenta associated with the fermionic variables are up to a constant 
factor identical to these variables themselves means that ip+ and i/>_ cannot be expressed 
uniquely as a function of configuration and momentum variables separately. It is therefore 
sensible to view ( |17D as a constraint and replace all occurences of 7r + and 7r_ accordingly. 
Hamiltons equations can then be written as 

dH dH dH . dH 

dp ox oip+ oip- 

so that the Poisson brackets for even quantities y and odd quantities 9 have to be: 

r -i (dyx dy 2 dy 2 dy x . dy x dy 2 . dy x dy 2 \ 
{1/1,1/2} = I "5 o t^— -7T- + l ^r, — h*- 



,„ , ( 89 dy dy 39 89 dy 89 dy 



dxidpi dxidpi dip +i dip +i dip_ i dip_ i 

m n^ (99 1 d9 2 89 2 d9 1 . d9 x d9 2 . d9 x 89 2 
{9i,9 2 \ = I 7: q h 7: q Zth — ~a~i Z: 



dxidpi dxidpi dij) +i dil) +i dt/>_ . dt/>_ . / 

where summation over z is implied. The result is the same as in if the constraint flTTD is 

used to replace the occurences of 7r + and 7r_ in the formulae given there. 

We can then calculate the brackets between the dynamical variables and find the only non-zero 

elements: 

{xi,pj} = dij] i{i/> +i , if> +j } = Sij] = Sij. 

From these we can finally determine the following algebra relations between all the conserved 
quantities: 

{H, H} = {H, Q} = {H, Q} = {H, R} = {H, L} = 0, 

i{Q, Q} = i{Q, Q} = 2H, {Q, Q} = 0, {Q, R} = Q, {Q, R} = -Q, 

{L, Q} = {L, Q} = {L, R} = 0, {L u Lj} = e ijk L k , {R, R} = 0. 



5 The Coulomb problem 

The supersymmetric Coulomb problem, characterized by an l/\x\ repulsive bosonic potential, 
can now be obtained from ([I]) by setting 

U(x) = — £, (18) 
\x\ 2 

or equivalently by choosing W(|sc|) = 2a | ai| ^ in (|12|)-(|l4|). It follows that the equations of 
motion now read: 

x 3 25 / 7 \ 

X = -r-r; -i-r—r (^+(^_-x) + (^ + -x)V»_+(^+-^_)«--(^+-k)(^_-Sb)sb) (19) 

JC 2 x 5 \ 2 / 
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V>+ = ~^U--h^--x)x) (20) 
\x\2 V * / 

= (21) 
| a? 1 2 \ 2 / 

We will proceed to solve these equations for a Grassmann algebra with two generators in the 
next section. 

The formulae for all the conserved quantities need not be written out again for our special 
choice of U but it is worth mentioning the explicit form of the Hamiltonian: 

1 1 2^ f 3 \ 

H = -x x + — ( V+-0- - -(t/>+-£)(V>--£) ) • (22) 

^ \ x \ \X 2 V ^ / 

In the purely bosonic case the Coulomb potential admits a further well-known conserved quan- 
tity, the Runge-Lenz vector, given by: 

V RL = [L, x] - x 

Unfortunately this extra conserved quantity does not seem to exist for the supersymmetric 
case. We will outline a proof for this in the appendix. 



6 Explicit solutions 

We will now solve equations (p^9|)-(pl|) for the case of a Grassmann algebra with two generators. 
In other words, we assume that all physical quantities take their values in an algebra that is 
spanned by two Grassmann-odd elements £i and £ 2 > satisfying the relations 

£? = & = o, 66 = -66- 

Using a method similar to the one deployed in || we split the bosonic variable x into components 
according to the number of generators involved: 

x(t) = x°(t) + x ij (t)i^j = x°(t) + x\t). 

Taylor-expanding the potential function U then gives us 

U(x) = U(x°) + VU(x°)x\ 

Notice that all quantities are vector-valued. We will indicate components, if necessary, by lower 
indices. 

The fermionic variables are time-dependent multiples of 6 and 6 since for the particular 
Grassmann algebra chosen there can be no product of three or more generators. Therefore in 
products of bosonic and fermionic quantities only bosonic terms of zeroth order in the generators 
contribute. Writing 

tf_(t) = #i(t)6 + i£(*)6 
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both fermionic variables decompose into two terms with real-valued time-dependent coefficients 
iM and respectively, where I = 1,2. 

The equations of motion for the coefficients are similar to (|20|) and fl2T|) — the only differences 
being that if) + and ?/>- have to be replaced by the component functions i/^ and if)[_ and x has 
to be replaced by x°. 

While the fermionic equations thus do not look much different after the decomposition the 
bosonic equation (|T9| ) splits in two, namely 



*° = 1% ( 23 ) 
\x u \^ 

Z 1 = Xl ~%f )& ° +a(x ^ + ,^_), (24) 

where 

a(x°^ + ,^_) = (V'+(V'-^ ) + (V' + ^ )V'-+(V' + -V'-)^ -^(^+^ )(V'-^ )* 

does not contain x 1 . 

Following || we will now systematically solve these equations by working "upwards" from the 
bottom layer of the system (^) through the fermion equations (|20D and fl21~D to the top layer 

m- 

6.1 The bottom layer bosonic equation 



The lowest order equation ( P5| ) is of course just the familiar equation of the Coulomb problem. 
Its solutions are the well-known hyperbolae in the plane orthogonal to the angular momentum 
vector L (if we set apart for a moment the special case L — 0). In order to make use of the 
solutions for the remaining equations it will be necessary to write them down in an explicit 
time-dependent fashion. To do this we utilize a method first devised by Moser in JTT| for the 



Kepler problem with a negative total energy and then extended to the positive energy case by 
Belbruno @. 



We can summarize these results adapted to our problem by saying that there is a diffeomorphism 
between the constant energy surface H = E° in the phase space of the Coulomb problem and the 
tangent bundle of the upper (or lower) sheet of the hyperboloid 7i specified in four-dimensional 
Euclidean space by Xq — X\ — X\ — X\ = 1. This diffeomorphism takes geodesies on 7i into 
solutions to ( p3|) requiring only a change of the time variable. 

More explicitly, if we parametrize a geodesic on TL by Xi(s); i = 0,. ..,3 and denote by 
Pi(s); i = 0, . . . , 3 its tangent, then 

x°(s) = P(s)(l+X (s))-X(s)P (s) (25) 
1 + X (s) 

transforms this geodesic into a solution of the Kepler problem for the energy E° = |, when we 
assume the following transcendental relationship between the geodesic "time" -parameter s and 
physical time t: 

t= f S Ix^s'^ds'. (27) 
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Solutions for arbitrary energy can then be obtained by the scaling 



x° = ^—x°, p° = (2£°)V, t = — ?— s-t. (28) 
2E° ' y v ; y ' (2£°)§ 

Every geodesic on W can be mapped using an appropriate SO(3, 1) -transformation into 

X(s) = (cosh s cosh /3, sinh s, 0, — cosh s sinh/?). 

For the Coulomb problem this equates to choosing the plane of motion and the orientation of 
the hyperbola in it. (Note that for the specific parametrization chosen we have also fixed the 
origin of time.) Specifically we get 

a: = ( e + cosh s, -(e 2 -l)^ sinh s, 0) (29) 
= sinh, («?-!)* cosh, 
1 + e cosh s 1 + e cosh s 

where e = cosh /3 represents the eccentricity of the hyperbola. The relationship between s and 
t is given by 

t = e sinh s + s. (31) 
This equation cannot be inverted analytically. 

One can read off from ( p9"D that the motion is indeed hyperbolic, taking place in the Xi-X2-plane 
with the focal points on the xj-axis and the origin of time chosen such that t = at the point of 
closest approach to the potential centre. We just mention that L is oriented in (— ^-direction 
and (for constant energy) \L\ = (e 2 — 1)5; the scattering angle 9 is determined by esin| = 1. 
There are two special cases worth mentioning, namely e = 1 when the scattering angle 9 is tt 
and the angular momentum is zero and e — ► oo when 9 = and the equations describe free 
motion at an infinite distance to the center of the potential. 



6.2 The fermionic equations 

To solve equations (|20|) and fl2T|) for x = x° we can first go back to our general solutions (|10|) 
and flTj]) . They state in this circumstance that and ?/>- are linear combinations of x° and 
U(x°) with Grassmann-odd coefficients. (Note that we can safely assume E° > due to the 
repulsive nature of the Coulomb potential.) 

To analyse the fermionic movement it is therefore sensible to look at the motion of x° and U(x°). 
The first interesting aspect is that if we combine velocity and potential in a six-dimensional 
vector Z = (x°, U) then the motion of this vector takes place on S 3 , since 

x°-x° + U-U = 2E° = 1 

for our particular solution, but Xg = C/3 = trivially. This is in fact very similar to the one- 
dimensional problem where motion takes place on S 1 instead 0. 

If we project the motion into the x\— i^-plane then this two-dimensional vector describes a circle 
with center at (0, — e/(e 2 — 1)2) (see Fig. 1). The motion starts for t = —00 at the intersection 
of this circle with S 1 , running initially towards the origin but then bending towards the x%- 
axis and reaching its nearest point to the origin at (0, — (e 2 — l)V(e + !))• Then it turns back 
outwards and runs symmetrically to the other intersection point, reaching it at t = 00. 
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Figure 1: Projection of the bosonic motion into the x^—x^-plane (e = f) 



The projection into the U\ — [/2-plane shows a different picture. Due to the factor 1 / 1 cc 1 2 in 
(|18|) the motion starts and ends at the origin. It runs radially outwards along an asymptote 
determined by the angular momentum but then bends towards the £7i-axis, crossing it at time 
t = and returns to the origin in a symmetric way (see Fig. 2). 
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Figure 2: Projection of the bosonic motion into the U\ — U2-plane (e = |) 

When we come back to solutions ( |T0| ) and (|TTJ) we can immediately state that 

lim = lim = Qx° 
lim ifj_ = lim = Qx°, 

so that we know the asymptotic behaviour of these solutions. 

However, we can visualize the fermionic motion completely if we refer to the bosonic motion 
on S* 3 . For this it is more convenient to use the real-valued components and than the 
Grassmann-valued quantities. One advantage is that we can meaningfully define the length of 
these vectors: For example, 

s 1 = \(¥ + ? + \m 2 

is a conserved quantity arising from invariance of the component equations of motion under 
change of generators. Note that because we are dealing with real-valued quantities now, none 
of the two products is zero. 

To understand the close connection between the above-mentioned bosonic motion on S 3 and 
fermionic motion let us ignore the third components ifr+o and ?/>_„ for a moment, since, as 
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we will later see, their equations decouple. Combining the remaining fermionic terms into the 
vector-valued object 

we see that ty 1 lies in the plane spanned by the bosonic vector Z = (cc°, U) and one specific 
orthogonal vector Zj_ = (-U, x°) obtained by rotating Z by ninety degrees in the Xi — Em- 
planes for i—1,2. (We are neglecting the trivial third components here). The angle between 
\& 7 and the two bosonic vectors is entirely specified by the two supercharges Q and Q. (To 
be more precise, it is specified by their real-valued components Q 1 and Q 1 with respect to the 
two Grassmann generators. For simplicity, we will not indicate this difference in our notation 
assuming that the precise meaning is always clear from the context.) In other words, as the 
bosonic vector Z moves on S 3 , determined by the bosonic equation of motion, the fermionic 
vector simply follows this movement rigidly, its direction being fixed with respect to the bosonic 
vector and its orthogonal complement by the two supercharges. 

This result is very similar to the one-dimensional case where one finds similarly constructed 
bosonic and fermionic vectors corotating on S 1 . Still, there is a major difference in more than 
one dimension. This lies in the fact that Q and Q do not specify the fermionic solutions uniquely. 
The solution in one dimension can be derived by inverting equations (0) and (H), yet in two or 
more dimensions this inversion is no longer possible. Instead these two equations only specify 
the two angles between ty 1 and the two bosonic vectors Z, Z± but not more. This means that 
the general ty 1 (again neglecting third components) moves simultaneously on two regular four- 
dimensional cones with center at the origin and their symmetry axes Z and Z±, respectively. 
As Z and Z± move with time ty 1 follows rigidly attached to the surface of both cones and thus 
lies on the intersection of them. In the particular solution that we gave above this intersection 
corresponds just to one point (taken at fixed distance from the origin) - thus specifying this 
particular solution completely. However, in general the intersection of the two cones (and 5" 3 ) 
will yield a one-dimensional sphere (just as the intersection of two three-dimensional cones 
and S 2 is equivalent to a zero- dimensional sphere, i.e. two discrete points). All points on this 
sphere have the same combination of supercharges Q and Q, which means that there exists a 
one-dimensional space of degrees of freedom for the fermionic motion that is not determined 
by the supercharges. 

In fact, we will now derive a second solution to (p0|) and (pl|) that fixes the motion of ^l/ 7 on 
this one-dimensional space. Since the particular solution involving the supercharges described 
above lies entirely in the Z-Z^-plane we will now look for a solution in the plane orthogonal 
to it. (Remember, that the space we are considering is four-dimensional.) This plane can be 
spanned by the two vectors: 

Y = {xl-xl-U 2 ,U x ) 
Y ± = {U 2 ,-U x ,xl,-x\). 

A straightforward ansatz for the second solution is ^> = PY + PY±, where P and P are 
Grassmann-odd constants. Yet, this ansatz turns out to be too simple and does not yield a 
solution. Instead we have to allow for time-dependent coefficients of Y and Y±. In other words 
we make the assumption that 

V = \x{t)Y + \ 2 {t)Y ± , (32) 
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where the A, are Grassmann-odd functions of time which still have to be determined. Inserting 
this ansatz into (^) and fl2"T| ) we find that the equations of motion are satisfied as long as the 
coefficient functions \i(t) obey the condition: 

^ = n , 3 e ij -Aj • (33) 

2a |ar| 2 

This first order differential equation for Aj can be solved by: 

A x (t) = Pcosuj(t) - Psinw(t) (34) 
A 2 (i) = Psincu(t) +Pcoscu(t), (35) 

where P and P are two Grassmann-odd constants and 

/•* 1 1 /•«(*) 1 

o;(t) = / — = — / r ds. (36) 

Jo 22 2 22 jo (l + ecoshs)2 

This integral cannot be evaluated analytically but it can be written in terms of the standard 
first elliptic integral F(s, k) as 



„ M /x / 2e V 



We can immediately read off from ([3(J) that u;(i) is a monotonically growing function that 
converges to a constant as t — > ±oo and that has only one zero at t = 0. The asymptotic 
behaviour of uo can be described by the formula 

7T 

Uoo = lim u;(£) = — lim u)(t) = — rP - 
t^+oo v ' t-»-oo v 7 2e5 _! 

where Pi{k) is the Legendre function of first kind. 

There are two interesting limiting cases, namely e = 1 and e — > oo. Because P_ 1/2(1) = 1 and 
P_!/ 2 (0) is a finite constant we find that u — > | in the first and — > in the second case. In 
other words, uj^ is a monotonically falling function of e, taking values in the interval (0, |]. 

Going back to equations (|34D and (^) and inserting them into our ansatz fl32l) we find thus the 
result: 

* = P (cos lu (t)Y + sin to (t)Y ± ) + P(- sin cu(t)Y + cos to (t)Y ± ) (37) 

which can be interpreted as follows (substituting the real components ty 1 , P 1 , P 1 for the 
Grassmann- valued quantities, if necessary): In addition to motion in the Z — Zj^-plane the 
fermionic vector can also move in the plane orthogonal to it. As in the previous case it is 
rigidly attached to and corotating with two orthogonal bosonic vectors. The motion of these 
two vectors, however, is slightly more complicated than before. They are not rigidly connected 
to Z and Zj_ in a vierbein but rather have a time-dependent phase (unless the eccentricity 
tends to infinity so that our formulae describe a freely moving particle at infinite distance to 
the origin, in which case u(t) = 0.) 

This phase changes in a continuous fashion between — c^oo and u^, where depends on the 
eccentricity of the bosonic solution and takes its maximum value | in case e = 1, i.e. in the 
absence of angular momentum. 
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Thus the second solution describes a circular motion orthogonal to the first but with a time- 
dependent phase difference with respect to the rigid frame provided by the bosonic vector Z 
and its three chosen orthogonals. We can think of this motion as parametrizing the motion on 
S\ that we mentioned above. 

So far we have neglected the third components of and those in the chosen direction of 
angular momentum, for which the bosonic movement is trivial (i.e. X3 = U3 = 0). It was one 
of the results from || that triviality of the bosonic solution does not extend to the fermionic 
solution and the same can be observed in the multidimensional case. 

Since the bosonic motion takes place in the Xi-a^-plane and therefore 23 = equations (^H) 
and (PH) simplify considerably: 

2^-25 

^+3 = ip-s = —^i^+s- (38) 

|a3 u p |a? u |2 

Incidentally, up to a factor of two these are the same equations as fl33|) and therefore have 
(again up to a factor of two) the same solutions: 

iP +3 (t) = Ocos2u(t)-dsin2u(t) (39) 
if)_ 3 (t) = Osin2uj{t) + Ocos2u{t), (40) 

where O and O are two further Grassmann-odd constants and u(t) is the same function as 
before. The only difference consists of the extra factor of two and is responsible for a slight 
change in the asymptotic behaviour: Whereas the maximum phase obtained for t — > ±00 
still converges to zero for e — > 00 implying that = O and ips = O are both constants, 
this maximum phase now tends to tt rather than ~ for e = 1 meaning that for zero angular 
momentum and describe a full circle. 

So although there is no bosonic motion in the third direction (to lowest order) the fermions 
somehow "see" the bosonic potential and behave accordingly. In fact, we can think of ([38D as 
fermionic equation of motion for a one-dimensional Coulomb potential - the only difference is 
that the strength of this potential is quadrupled with respect to the original problem. To see 
this we have to interpret the right-hand sides of (J53) as TU'(x)tp T where U'(x) is given by 

U'(x) = % U(x) = -K- 
X2 X2 

Thus the imagined one-dimensional bosonic potential has to be 

V(x) = -U 2 = - = AV. 
2 x 

This can be seen as an explanation for the extra factor of two appearing in the solutions ([39]) 
and 

We want to finish this section by a short discussion of the newly found fermionic constants P, 
P, O and O since there is a major difference between them and the supercharges Q and Q. 
Formally we can write all four constants as conserved quantities by simply inverting equations 
(|32|), (|34|), (|35|), ( |39"D and fliPp . However, they cannot be seen as an original symmetry of 
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the Lagrangian, in the sense that they cannot be derived as Noether charges by a variation 
just involving the dynamical variables. On the contrary, to derive those charges from the 
Lagrangian we have to presuppose a certain knowledge of the system, e.g. the orientation of 
angular momentum, to find the correct variation. As an example, O and O can be formally 
found as conserved quantities using the following variation: 

5i[) +s = — e cos 2u(t), 5ip-3 = — esin2co>(t) 
5ip + 3 = esin2u;(t), Sips = — ecos2u;(i). 

Evidently, neither function on the right hand side is a function of the dynamical variables only, 
instead we find a rather complicated function of time involving u(t). In addition this variation 
would have to be different if we would alter the initial data of the bosonic motion, for example, 
by making a different choice for the orientiation of angular momentum. Thus P, P, O and O 
should be seen as "lesser" charges, not comparable with the two supercharges Q and Q which 
are genuine Noether charges as we have seen above. 

There are, however, some interesting observations about the interrelation of these extra 
fermionic quantities to be made. Firstly, one can read off from (^9]) and (|40|) that the product 
S = iip+3ip-3 or, more generally, S = i(L° -ip + )(L is conserved, where L° denotes the 
bottom layer angular momentum. This is a consequence of the symmetry 

5if> + = e{L°-if>_)L , 5<4>- = e(L°-il> + )L . 

Since we have restricted ourselves to two Grassmann generators we can replace L simply by 
L and so find a genuine symmetry of the Lagrangian, involving only the dynamical quantities. 
Moreover, between all six fermionic constants mentioned there then exists the remarkable re- 
lation 

i(QQ + PP + 00) = 2ER, 

where, as we have seen, the bosonic constants R and E are genuine Noether charges for arbitrary 
Grassmann algebra. 

6.3 The top layer bosonic equation 

We finally have to solve equation ( f2~3| ) for the bosonic top layer, concentrating first on the 
homogeneous part of this equation, i.e. ignoring a(x°, if) + , ?/'_). 

In the one-dimensional case there were two solutions to the analogue of this equation, namely 

ft 1 

x 1 = cx°; x 1 = cx° / —p—dt'. 

J to (x u r 

Whereas the first term is still a solution to the three-dimensional problem, the second one is 
not. A more fruitful generalization from the one- to the multidimensional case is given by the 
idea that solutions to the top layer equation describe variations of the bottom layer system 
with the free constants of that motion. We will explain this idea in more detail below. 

To find an explicit form for the solutions we try the following ansatz 

x\t) = f 1 (t)x (t) + f 2 (t)x (t). (41) 
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Inserting this ansatz into (j2"4|) and using the bottom layer equation (|23"D we can evaluate both 
sides and write them as time-dependent linear combinations of x° and x°. Comparison of the 
coefficients then yields the two relations 

m+j^m = (42) 

2fi(t) + f 2 (t) = 0. (43) 
Integrating equation ([43]) and inserting the result 

/ 2 (f) = -2/i(f) + c (44) 
into fl42| ) we get the following equation for fi(t): 

h(t) = j x \^(f l (t)-2c). (45) 

We can recognize that the homogenous part of this equation is nothing else than the bottom 
layer equation and this means that we can write its two solutions as the two component solutions 
of ( p3|) given in (p9|). Since the inhomogenous solution is just the constant 2c we can write 



fi(t) = Cix1(t) + C 2 x° 2 {t) + 2c = d(e + coshs(t)) - C 2 (e 2 -1)3 sinhs(t) + 2c. 
Applying this result to (f44"l) we find that 

f 2 (t) = -2C x X x {t) - 2C 2 X 2 (t) - 3d + d, 

where 

X x (t) = fox1(t')dt' = \es{t) + lesinhs(t)coshs(t) + (e 2 - l)sinhs(t) 
X 2 (t) = x° 2 (t')dt' = -(e 2 -l)^(coshs(t) + iesinh 2 s(t)-l). 



Inserting these results back into our ansatz (|41|) we find that 

aj 1 ^) = Ci (x^x^t) - 2X 1 (t)x°(t)) + C 2 (x° 2 (t)x°(t) - 2X 2 (t)x°(t)) (46) 
+c (2x°(t) - 3tx°(t)) + cx°(t). 



There are still two solutions to (^) missing. As all solutions found so far lie in the plane of 
motion spanned by x°(t) and x°(t) we have to look now at solutions in the third direction, i.e. 
for xl(t). The equation of motion simplifies here to 

= pp^W 

since x®(t) = 0. Again this is just the bottom layer equation of motion, solved by x±(t) and 
x 2 (t). The general solution to the homogenous part of ( p4|) is thus given by (^6|) and 



x 1 ^) = C 1 x° 1 (t)L° + C 2 x° 2 (t)L°. (47) 
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When we look for an interpretation of these solutions the analogy with the one-dimensional 
case will be - as mentioned - quite fruitful. In one dimension the solution to the homogenous 
part of the x 1 -equation could be written as a sum of variations: 

ll(i ) = Cl g (( ) +C2 ^ (i) . 

The same turns out to be true in the three-dimensional case. 

To start with, the term with coefficient c in fl46|) can be clearly interpreted as variation of our 
bottom layer solution with initial time to : 

dt 

Variation with energy can also be found in the solution for x 1 as can be seen in the following 
way: While solution Q2§D is valid only for energy E = | we can obtain a solution for arbitrary 
energy E by scaling both x° and t according to (p8[). The variation of x° with energy is then 
obtained by comparing x° E (t) with x° E+dE (t): 



6x° E _ dx° E + dx% dt 



5E dE dt dE 2E 2 
Evaluating this result at E — ~ yields 



x° + 3(2E)h E x° E . 



5x 



5E 



-2x° + 3tx l 



leaving us with the term in (jig ) with coefficient c. It is quite important to realize that we 
have to compare two solutions of different energy at the same physical time t, i.e. we must not 
change clocks and use t E as a parameter for the solution x° E+dE . 

These two variations comprised the homogenous part of the solution to the one-dimensional 
problem. For the Coulomb problem there are four more independent ways to vary the bottom 
layer solution. One of these consists of rotating the hyperbolic orbit in the plane of motion by 
a small angle <j) around the angular momentum axis. The variation is thus given by 

5X Q Q 

"^0" = ( — X 2i X li 0)- 

This variation coincides with a linear combination of two terms in (^6|), namely 
Sx° e 



e 2 - 1 



We can thus interpret the term with coefficient C2 in (46) as a combination of rotation in the 



plane of motion and variation of the initial time parameter. 

Another way to vary x° is to change the eccentricity of the hyperbola which is explicitly included 
as a parameter in the solution (p9| ). We find 

Pt^Q Pic / ^,0^,0 ^o„o \ 



Sx° dx° dx° ds 



2 1 11 



— i— - — _ = 1 - z x 11 

5e de + ds de { L° 3 ' L° ' 
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where we have used that 



ds I de \ sinh s 

de \ds I 1 + ecoshs 

t=const 

We again find this variation as a linear combination of two terms in fl46|) , namely: 

5x° 1 



x\x Q - 2Xi x°) - e (2a; - 3tx 



be e 2 — 1 

and so the the Ci-term in ([5]) describes a combination of variation with eccentricity and energy. 

It is now to be expected that the two remaining terms of the homogenous a^-solution are 
connected to a tilt of the plane of motion or, equivalently, to a change in the orientation of 
angular momentum L°: 

5L° = 5L° iei + 5L° 2 e 2 , 
here e% and e 2 represent unit vectors in the respective directions. We find that 

5x° 5x° 



x i e 3> TTo = X 2 e 3 



6L1 1 d ' 5L° 2 



which are evidently proportional to the two terms in (47) since L° is oriented in the third 
direction in our particular solution to the bottom layer system. 

Thus we have now derived and found an interpretation for all the terms of the homogeneous 
solution to the top layer system. As was the case for one- dimensional systems this solution 
describes all the possible variations of the bottom layer with all possible parameters describing 
this layer. In our case we need six parameters to describe the motion of the original Coulomb 
problem (initial time, energy, eccentricity, orientation of major axis and two parameters de- 
scribing the plane of motion) and correspondingly there are six independent solutions to the 
top layer. 



The last step is now to find one particular solution to the inhomogenous equation fl24}) . The 
general idea is to take linear combinations of the inhomogeneous terms a.j, integrate them over 
time and multiply them by one of the solutions x\ om of the homogeneous equation, then sum 
over all solutions: 



-y /"* III 

X inhom = X hom.a / fa,i(t )dj(t )dt . 

1 Jo 



a=l 

We determine the coefficient functions f a ^ by two conditions: Firstly, in taking the time deriva- 
tive we want to get 



t 1 

inhom 



~, Jo 



a=l 

i.e. all the terms involving derivatives of the integrals must cancel out each other. Then for 
the second time derivative we arrive at 



inhom 



n t / n \ 



a=l " u \a=l 
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Because the x\ om a -terms satisfy the homogeneous equation the first sum will give the homo- 
geneous part of (0). For the next term to give the inhomogeneous part we need the second 
condition that the sum in parenthesis is equal to the unit vector e\ 

We will now demonstrate this procedure: For simplicity we begin with the third component of 
x 1 . Since it decouples from the other components we can safely assume f a ,i — fa, 2 — 0> 
writing f a instead of f a 3 we can make the ansatz 



,1 

J inhom,3 



, =x° 1 / h(t')a z (t')dt' + x° 2 / f 2 (t')a 3 (t')dt'. (4f 



Our first condition then yields Xi(t) fx(t)a 3 (t) + x^it) f 2 (t)a 3 {t) = and thus we can set fi = cx\ 
and f 2 = —cx\. Then we find 

x] nhomt3 = c (x? jT x°(t')a 3 (t')^ " ^ *S(f MO^) + <44 ~ 44H(t) 
= j^-c (xl f x°(t')a 3 (t')dt' - x° 2 f x\{t')a 3 {t')dA - cL° 3 a 3 (t) 

X \ J •/ / 

and so by equating c = -4 we end up with the required result, namely 

^3 

%inhom,3 = |^0j3 X 3 + a 3(/)- 

Since the equations for the first and second components are coupled we need all four remaining 
homogenous solutions to give the correct inhomogeneous term. We just write down the result: 



inhom " jjo 



*° f Q ax (^44 + 2§ 3 (L° 3 ~ ±\xl) - 3tx°\ (49) 

+a 2 (2^x° l x° 1 - 2^(L° 3 + x°x°) + 3tx° 1 ) dt' 
\ L 3 L 3 J 

+ (2x° - 3tx°) (-a x x% + a 2 xf) dt' 



0^,0 ™0™0 

Ju Ci Ju 



+ (x lX ° - 2X lX °) j o L^-a 2 (l + 



2^1 ■ 



TO ) 

3 ^3 , 

0^0 ™0„0' 
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+ (x° 2 x° - 2X 2 x°) J L(l - ^) + a 2 



3 ^3 



dt' 

dt' 



Because we know both the bottom layer bosonic and the fermionic quantities contained in a 
we can now in principle insert these into the integrals and evaluate the solutions as explicit 
functions of time. However, since the knowledge of these functions does not provide new 
insights, we refrain from writing them down here. 

In summary, the complete solution to the top layer bosonic equation is given by the sum of the 
homogeneous solutions ([46]) and (f|7|) and the inhomogeneous solution given by ( flgf) and d49|) . 
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7 Discussion 



Our aim in this paper has been not only to analyze and solve a fascinating problem in supersym- 
metric classical mechanics but also to stress the relations between supersymmetric mechanical 
models in one and more than one dimensions. An example is that the fermionic variables should 
really be thought of as components of a 2n-dimensional vector \1/ which moves on a (2n — 1)- 
sphere. In the one- dimensional case this means that all fermionic vectors move on circles - 
which is consistent with our findings in ||. Furthermore, the motion is coupled to the motion 
of a purely bosonic vector, made of the dynamical variables x and U, by the supercharges. 
While these completely determine the fermionic motion in one dimension we have seen that the 
situation in higher dimensions is more complicated and the supercharges only prescribe a geo- 
metric subspace for the true solution. Motion on this subspace is then no longer characterized 
by the rigid connection of bosonic and fermionic motion found in the one-dimensional case but 
exhibits a rather interesting time- dependent phase. 

We have shown that for rotationally invariant scalar potentials there will always be a conserved 
angular momentum so that to lowest order bosonic motion takes place in a plane. It is not 
unreasonable to conjecture that for every direction orthogonal to that plane we will have to 
solve a one-dimensional problem for the fermionic quantities corresponding to that direction 
- as was the case for the Coulomb problem. This leads us to another generalization from the 
one-dimensional case, namely that trivial solutions to the bosonic equations of motion do not 
imply that the corresponding fermionic equations can only be solved trivially. In fact, for our 
particular problem we found that the third components of the fermionic quantities still "felt" 
the presence of the bosonic motion in the plane orthogonal to them. 

One of the most fruitful ideas for the higher bosonic layers which can be applied in any di- 
mension turned out to be that solutions for these layers are variations of the solutions of the 
classical problem. While in one dimension every solution can be conveniently described by the 
energy and some choice for the initial time, in higher dimensions more of these parameters are 
needed - in case of the Coulomb problem we have conveniently related them to the elements 
describing the classical orbit. So the view taken in || for the one-dimensional case, namely that 
supersymmetric dynamics captures information over a whole range of energies of the system, 
can be confirmed in a more general sense for our higher dimensional model: The bosonic (but 
Grassmann-valued) part of the solution describes not only one particular classical solution (by 
its component of order zero) but also includes all possible variations of it with all possible free 
parameters or, in other words, a bosonic solution corresponds not just to a point in the real 
parameter space but to something like a fuzzy subset. 

For reasons of simplicity we have described in this paper solutions for a Grassmann algebra 
with two generators. From our experience with the one-dimensional model and the clear links 
that it exhibits to the three-dimensional system we have studied in this paper, we conjecture 
that the picture for a larger number of generators is similar to that in the one-dimensional case. 
There we were able to show that solutions to the fermionic equations corresponding to three 
generators consist of two parts: One that looks like the first order solutions themselves (i.e. 
movement on a circle in one dimension) and one that consists of the variations of these solutions 
with energy and initial time. Similarly, the bosonic equations for four generators are solved by 
the first and second variations of the lowest order bosonic solutions and first variations of the 
interaction term. 
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It is therefore not unreasonable to assume that the general picture in three dimensions is not 
much different and that the higher order equations of motion can be solved by a suitable 
combination of (first and higher order) variation terms with respect to the free parameters of 
the physical model in question. 



Appendix: Runge-Lenz- Vector 

In this section we want to outline a proof that the supersymmetric Coulomb problem discussed 
in this paper does not admit a conserved Runge-Lenz vector. In the classical problem this 
conserved quantity cannot be constructed as a Noether charge but arises as a consequence of 
the hidden 5*0(3, l)-symmetry of the problem, explicitly recognizable in our construction of 
the bottom layer solution to the equations of motion in section |6.1| . This symmetry must be 
broken in the supersymmetric case. 

We start with the assumption that there is a conserved Runge-Lenz vector in the supersym- 
metric problem. Setting all fermionic variables to zero (for our choice of Grassmann algebra 
this can be accomplished by choosing Q = Q = P = P = O = O = 0) we expect to get the classical 
result for this vector, namely 

Vrl = [L,p]-x, (50) 

where L is given just by its bosonic part [a;,p]. (We will use p = x for this section). 
Notice that the Hamiltonian and the angular momentum vector also have this property. Of 
course, here all quantities are Grassmann-even rather than real, but this makes no difference 
and, of course, fl50D is a conserved quantity if we set all fermionic terms in the Lagrangian ([!]) 
to zero. 

Returning now to the general case with non-trivial fermionic variables we take the time- 
derivative of (^) and find (still taking the classical result for L): 



d 3 . 25 



-x (2(ip + -p) (V>-- x) + 2(il> + - x) (V?--p) + *!>-) (p-x) - 7 - {il>+- x) (V>- x) (p-x) 

The result is clearly non-zero but then we can expect the supersymmetric version of the Runge- 
Lenz vector to have an additional fermionic piece. However, the derivative of this piece should 
yield exactly the same result as above with an extra minus sign so that the derivative of the 
total will be zero and hence V rl a conserved quantity. We will now try to construct this 
fermionic piece. 

We start by writing our result for the bosonic part of V rl in components, so that the linearity 
in p, and can be seen in a more explicit way: 



d T rU 3.2a / .. , , 3 



_ y*. = — i — I [5 lj 5 kl — 5 kl x i x j +5 lk x j x l +5 ]k - i - 1 
dt 2 I a? 1 2 V 2 



X X 
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We now try to assemble all terms whose derivative could possibly contribute to this result. For- 
tunately, there are restrictions as to which terms to consider that shorten the list of expressions 
to look at. 

Firstly, we need only to analyze terms quadratic in the fermionic variables: Since Vrl is even 
their number has to be divisible by two, all terms containing no fermionic quantities at all can 
be assumed to be contained already in fl50|) and if there were four or more fermionic variables in 
Vrl the time derivative would still contain the same number of these variables in contradiction 
to flSlf). Furthermore since our system is symmetric under the exchange ij)^. — > —if)-, tp- —>■ tp+ 
we have to include only terms subject to this symmetry. This leaves us with multiples of the 
terms 

A) iipiipi. 

B) itp l + ip J + + iip l _ip J _. 

When we now classify all terms according to the type of their fermionic part and the even or 
odd number of p'-variables involved in each one of them, then the time derivative can be seen 
as acting on these classes as follows: 

type B, p l even 
type B, p l odd 
type B, p l odd 
type B, p l even. 

Here p l even/odd means the number of p'-variables contained in a term is even or odd, respec- 
tively. 

Since we want to derive a term of type A with an odd number of p h s, namely just one, we can 
read off from the summary above that we need not consider combinations of type II or III, since 
their time derivatives do not yield such a term nor can they compensate the B-terms that will 
arise from combinations I or IV. This means we only have to consider these latter combinations. 
However, even those are not unproblematic: As mentioned, they give rise to fermionic terms 
of type B which we do not want in our derivative, so we must hope that we can cancel these 
terms against each other. 

In the next step we make - preliminarily - the assumption that the number of p'-variables 
involved is either zero or one. We will later see that allowing for higher powers of p l will not 
change anything. 

Up to this point the allowed combinations left are 

a)^+^-, b) (i-0+V J + + iiptip-)p k ■ 

We now have to contract these with the remaining natural tensors, namely 5^ and x l , to 
leave only one free index. For terms of type a) we have the following options: x % x^x k , x l S^ k , 
x j S lk , x k 5 lj . 

For terms of type b) we have some extra restrictions: First, since they are antisymmetric under 
the exchange i <-> j we need not to consider -terms that contract the two fermionic variables. 
Second, if we include any x-term the time-derivative will automatically yield an additional p l - 
term so that the result would be quadratic in the p s - a case which we have excluded above. 
Therefore the only option left for combination b) is S^ k 5 d . 



dt 

I) type A, p l even — ; 

II) type A, p l odd — ; 

III) type B, p l even — ; 

IV) type B, p l odd — ; 



type A, p l odd + 

type A, p l even + 

type A, p l even + 

type A, p l odd + 
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We still have to determine the prefactor of each term. It can be derived by dimensional 
arguments. As can be read off from the Coulomb-Hamiltonian the dimension of energy is 



length . From this we can calculate the dimension of every other dynamical quantity: 



X 


P 


</>+ 


</>- 


v RL 


length 


length -172 


length 174 


length 174 
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Since the bosonic part of the Runge-Lenz vector is dimensionless we have to multiply each 
candidate term for the fermionic part by the appropriate factor of r = \x\ to make it dimen- 
sionless, too. We end up with the list in the first column of the following two tables. These 
tables describe the action of the time-derivative on each candidate term. In the top rows we 
have specified all tensorial combinations that can possibly be generated by applying the time- 
derivative to the candidate terms in the first column. The table itself then specifies the correct 
numerical prefactors of each tensor in the time-derivative of a particular candidate term. For 
example, one can read off from the table entries for the third term that 



l-8 ik 8 jl ^ 



8 lk x j x l ) r zpHip+tp 3 



For comparison we have added the bosonic part of the Runge-Lenz vector and the tensor 
decomposition of its time-derivative in the last row. The result is the following: 
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We now have to choose appropriate coefficients for the candidate terms such that their sum 
equals the bosonic part of V rl (with an extra minus sign). Looking at the coefficients of 
fiiifijk anc j fitkfiji we can immediately see that terms 2) and 3) cannot occur since these tensorial 
combinations are not part of the bosonic Runge-Lenz vector and cannot be compensated by 
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any other term. Looking then at the second table we can conclude that term (5) cannot occur 
either - since it contributes a non-zero factor to both columns whereas V rl does not (and the 
terms 2) and 3) have to be zero as deduced above). 

The column for 5 l ^8 kl tells us that the coefficient of the fourth term has to be |v^2 - in 
contradiction to the column for 5^x k x l which says that the same coefficient has to be v2. 
Similar contradictions arise for our first candidate term when we compare the columns for 
5' ll x^x k , 6^ l x l x k , 5 kl x l xi and x l x^x k x l . This leaves us no other option than to conclude that 
none of our candidate terms can be added to the bosonic part of the Runge-Lenz vector to 
achieve a zero derivative. 

One loophole remains. We have excluded further terms of type b) on the grounds that they 
yield terms quadratic in momentum. We might hope, however, that in the right combination 
these quadratic terms can cancel each other. There are two possibilities, which can be found in 
the first column of the following table; all other tensorial combinations either yield zero or are 
up to minus sign identical to the two choices below via the index exchange % <->• j. On taking 
the time derivative of both terms the crucial issue is whether the terms quadratic in p l can be 
cancelled against each other, so only quadratic components are shown here: 
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Evidently, the quadratic terms do not cancel each other. We conclude that no term at most 
linear in momentum exists that could be added to V RL,bos to give a conserved quantity in the 
super symmetric model. 

We can of course ask the question whether the fermionic part of the Runge-Lenz vector could 
be quadratic or of some higher power in p k . We will now show that the answer to this question 
is no. The reason for this can be understood by looking at a general term like 

AV k i~ k ™ l f( r )p kl . ..p km iip\.ipi, (52) 

where A^ kl '" kml (I is a free index) is a tensor constructed from two types of building blocks, 
S pq and x p , and f(r) is some well-determined power of r. Since the tensor is of type A we can 
infer that the number m of momentum variables has to be even. Furthermore, to keep (|52l) 
dimensionless f(r) has to be r^" 1-1 ). Clearly for every m > 2 /(r) is a non-trivial function 
of r. This in turn means that the time-derivative in acting on /(r) will produce an extra p k 
variable. 

There are only two ways to solve this problem: Either all the terms in the time derivative of 
(j52|) containing extra momentum variables cancel each other - or we have to make use of yet 
higher powers of p k in Vmjerm. to compensate. However, those very terms will by the same 
principle give rise to another extra power of p k in the derivative and so on. For this series 
to stop at some point we must hope that all terms in the derivative of A i: * 1 '" km with m + 1 
momentum variables cancel each other. We find that 

^_ fj^,jki...k m l r \{m-l)pki p^ilp^lpi = - (^ i 3 fc l-- A "^ r §( m ~ 1 )p fc l p^iip^ipj 

2 
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where the dots indicate further terms containing less than m + 1 momentum variables, in which 
we are not interested here. 

As mentioned above A 1 ^ 1 "' kml consists of a sum of tensors built from 5 pq and x p terms. While 
the time derivative does not act on the first type of term, it acts on the second: 



Thus we can write 



^x p = -(5 pl -x p x l )p l . 
at r 



dt r 



where A 1 ^ 1 " mkm¥1 is a new tensor derived from A l i x " m in the following way, which can be 
seen as an algebraic version of the Leibniz rule: We replace one x p tensor in A l ^ kl " ,kml with 
(5 pfcm+1 — x p x km+1 ), repeat this step for all other x p tensors contained in A % i kl - kml - replacing 
always one tensor at a time - and then sum up every contribution. Therefore 

^ijki...k m l r ~(m-l)pki p km iip'^ip ^ = {j[ijki...k m k m +il _|_ }i(^ n _\^J^jki...k m l^,krn+x\ 

So, if we want a cancellation of all terms containing m + 1 momentum variables, we need: 

j^ijk 1 ...k m k m +il _|_ }^^ m _^j^ijk 1 ...k rn l^,k m+1 _ g (53) 

It remains to show that this equation cannot be satisfied. Starting with 
we find 

£j kl ...k m k m+1 i = _( m+3 )^ £ i^i _ _ _ ^ijiifWi^ . . . x km x l +. . .+x i x j x kl . . . x km 5 lkm+1 . 

Inserting these results into fl53| ) we find the necessary condition — |(m+7) = 0, i.e. m = —7 
which is a clear contradiction, so a term like ( |54D cannot occur. Replacing two x k terms by a 
5-tensor we arrive at 

A\ jkl - kml = VVx^ . . . x km ~ l 5 kml (55) 
to name but one possible choice. The corresponding A is then given by 

^jk x ...k m k m+1 l = _( m+1 )£i£f£ki . . . x k m -i 6 k m l~k m+1 + _ _ _ 

where the dots indicate further terms which contain less than m + 1 £ fc -tensors. Equation (^) 
then yields — ~(m+3) = 0, which is again a contradiction. 

One can now see that as we subsequently replace two x fe -terms by one 5-tensor, equation (p>5|) 
gives that — |(m+7— 2p) = 0, where p is the total number of x fc -terms exchanged as compared 
to (|54]) . But this shows that m has to be an odd integer in contradiction to our demand that m 
be an even number. As a result no type A tensor can have the maximal number of momentum 
variables in the fermionic part of the Runge-Lenz vector. 

We refrain from repeating the argument for type B tensors since the result is the same: The 
time-derivative in acting on these tensors increases the number of p fc -variables by one and there 
is no way to cancel all these extra unwanted terms against each other. We can therefore safely 
conclude that even the introduction of higher orders of momentum does not change the overall 
result, namely that a supersymmetric version of the Runge-Lenz vector does not exist. 
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